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Abstract 

Ultracold  atoms  confined  in  optical  lattice  potentials  provide  an  ideal  platform  to  address 
some  of  the  most  outstanding  questions  in  quantum  many-body  physics,  that  are  difficult  or 
impossible  to  study  using  conventional  solid  state  systems.  During  this  DARPA  program, 
progress  in  several  important  directions  has  been  made.  We  have  developed  cooling  techniques 
suitable  for  a  lattice  system  and  cooled  fermionic  lithium  atoms  confined  in  a  three-dimensional 
(3D)  cubic  lattice  to  unprecedented  low  temperatures,  very  close  to  the  Neel  temperature.  We 
employed  Bragg  scattering  of  near-resonant  light  to  observe  the  antiferromagnetic  correlations. 
We  have  mapped  the  phase  diagram  over  the  entire  BEC-BCS  crossover  for  a  trapped  3D  gas 
with  imbalanced  numbers  of  spin-up  and  spin-down  atoms,  and  have  also  observed  the  crossover 
from  a  ID  to  a  3D  spin-imbalanced  Fermi  gas.  We  have  investigated  dynamics,  rethermalization 
and  out-of-equilibrium  phenomena  in  strongly  correlated  lattice  systems,  as  well  as  in  disordered 
lattice  systems.  These  studies  greatly  deepened  our  understanding  of  quantum  phases  of  matter 
in  lattice  potentials. 


During  this  funding  period,  significant  progresses  have  been  made  by  our  team,  several 
milestones  have  been  reached,  and  new  research  directions  have  been  identified.  The  following 
is  a  summary  of  our  research  achievements  made  possible  by  the  DARPA  funding. 

3D  Fermi-Hubbard  Model  and  antiferromagnetism 

The  Hulet  group,  in  collaboration  with  theory  groups  led  by  Huse,  Scalettar  and  Trivedi, 
have  used  Bragg  scattering  of  light  to  detect  antiferromagnetic  correlations  in  a  three- 
dimensional  (3D)  optical  lattice  with  densities  near  1  atom/site.  This  was  accomplished  by 
preparing  the  atoms  at  temperatures  of  1.4  TN,  where  TK  is  the  Neel  temperature.  This  is  a  factor 
of  two  below  the  previously  reported  lowest  temperatures  for  atoms  in  a  lattice.  By  comparing 


the  Bragg  signal  to  quantum  Monte  Carlo  (QMC)  calculations,  they  realized  thermometry  in  the 
lattice  that  is  accurate  at  the  10%  level.  The  closed  black  circles  in  Fig.  la  are  the  trap-averaged 
Bragg  signal  for  several  values  of  U/t,  where  U  is  the  onsite  interaction  and  t  is  the  nearest 
neighbor  tunneling  rate.  The  colored  lines  in  Fig.  la  are  the  results  of  QMC  for  several  values 
of  scaled  temperature  Tit.  Figure  lb  shows  the  calculated  Bragg  signal  vs.  Tit  as  open  diamonds, 
while  the  grey  area  is  the  uncertainty  range  for  the  measured  Bragg  signal.  The  data  fit  best  to  T/t 
=  0.50.  In  comparison,  T^/t  maximizes  at  0.36. 

The  development  of  a  novel  “compensated”  optical  lattice  was  essential  for  the  success  of 
this  experiment.  The  compensated  lattice  uses  blue-detuned  anti-confining  beams  superimposed 
on  each  of  the  infrared  lattice  beams  to  offset  the  inhomogeneous  confinement  envelope  of  the 
bare  lattice.  Compensation  provides:  1)  a  way  to  adjust  the  density;  2)  cooling  (or  mitigated 
heating)  in  the  lattice  by  evaporation;  and  3)  a  flattening  of  the  inhomogeneous  band  structure 
produced  by  the  lattice  confinement. 
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Figure  1.  a)  Closed  circles  are  the  trapped -average  Bragg  signal  Sq  (related  to  the  static  structure  factor)  for  Q  =  n, 
a  reciprocal  lattice  vector  for  antiferromagnetic  order  and  the  open  circles  correspond  to  an  arbitrary  Q,  which  is  not 
Bragg  enhanced.  The  colored  lines  are  the  result  of  QMC  for  the  various  temperatures  indicated,  b)  QMC  (open 
diamonds)  and  the  uncertainty  range  of  the  measured  Bragg  signal  showed  in  grey.  The  best  fit  to  the  data  is  T/t  = 


0.50. 


The  antiferromagnetic  phase  is  also  expected  to  be  a  Mott  insulator.  With  exactly  one 
atom  per  lattice  site,  the  compressibility  of  the  Mott  phase  goes  to  zero.  They  explored  the  Mott 
phase  using  in-situ  imaging  of  the  density  distributions.  An  example  of  the  column  density 
distribution  directly  recorded  in  these  images  is  shown  in  Fig.  2a,  while  the  inverse  Abel 
transform,  which  reconstructs  the  three-dimensional  distribution  is  shown  in  Fig.  2b.  These 
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distributions  have  the  characteristic  flat-top  of  the  Mott  phase.  A  Mott  plateau,  such  as  this  one, 
had  not  been  clearly  resolved  in  the  past  for  fermions. 


Since  this  data  was  obtained,  they  have  improved  the  stability  of  the  compensated  lattice.  Better 
stability  will  enable  us  to  better  optimize  the  compensated  lattice  for  cooling  to  lower 
temperatures.  If  successful,  they  will  explore  the  critical  regime  at  the  phase  transition,  and 
obtaining  even  lower  temperatures  may  open  up  the  possibility  of  observing  d- wave  superfluidity, 
which  has  been  conjectured  to  provide  the  mechanism  for  high  temperature  superconductors. 


Figure  2.  In-situ  images  of  a  3D  optical  lattice  for  various  interaction  strengths,  a)  column  densities;  and  b) 
reconstructed  three-dimensional  distributions  showing  Mott  plateaus  for  stronger  interactions.  The  large  number  of 
atoms  in  the  Mott  phase  is  a  result  of  the  flattening  of  the  lattice  band  using  the  compensated  lattice.  The  atoms  in 
the  Mott  phase  develop  antiferromagnetic  correlations  at  sufficiently  low  temperature. 


Spin-imbalanced  two-component  Fermi  gas 

Hulet' s  group  have  also  made  great  progress  on  the  experimental  study  of  spin- 
imbalanced  fermions.  They  have  mapped  the  phase  diagram  over  the  entire  BEC-BCS  crossover 
for  a  trapped  3D  gas  with  imbalanced  numbers  of  spin-up  and  spin-down  atoms.  A  manuscript 
reporting  these  results  is  nearing  completion.  They  have  also  now  observed  the  crossover  from  a 
ID  to  a  3D  spin-imbalanced  Fermi  gas.  This  experiment  uses  a  two-dimensional  optical  lattice 
to  create  a  bundle  of  ID  tubes.  As  the  lattice  depth  is  decreased,  the  coupling  between  tubes  is 
enhanced,  leading  to  a  crossover  from  a  ID  to  a  3D  phase  diagram.  They  have  found 
experimentally  that  the  location  of  the  dimensional  crossover  is  universal,  such  that  data  for 
different  interaction  strengths  and  lattice  depths  cross  from  ID  to  3D  at  the  same  value  of  f/sb, 
where  t  is  the  hopping  strength  between  tubes  and  Sb  is  the  binding  energy  of  the  pairs  in  the  ID 
tube.  The  elusive  FFLO  state  is  predicted  to  be  stabilized  near  this  crossover,  so  their  next 
experiment  will  be  to  search  for  the  expected  domain  boundaries  that  are  potentially  observable 
in  a  time-of-flight  experiment. 
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Several  theory  groups  have  worked  on  this  topic.  Pu's  group  considered  an  attractive 
two-component  Fermi  gas  with  Raman-induced  spin-orbit  coupling,  and  showed  that  under 
certain  situations,  the  Cooper  pairs  naturally  possess  finite  momentum  as  a  result  of  the  interplay 
between  spin-orbit  coupling,  spin  imbalance  induced  by  the  effective  Zeeman  field,  and 
interaction.  Recently,  Bolech's  group  has  extended  the  continuous  matrix  product  states  (cMPS) 
to  nonrelativistic  fermions  and  shown  that  the  cMPS,  correctly  describes  the  ground-state 
superfluid  and  magnetic  properties  of  interacting  Fermi  gases  in  ID,  as  those  realized  in  Hulet's 
lab.  This  includes  the  signatures  of  a  partially  polarized  superfluid  regime  (FFLO),  in  agreement 
with  the  large  amount  of  theoretical  and  experimental  work  from  the  DARPA-team  effort.  The 
new  type  of  ansatz  promises  to  be  ideally  posed  to  be  able  to  describe  atomic  gases  in  optical 
lattices  economically  but  without  making  a  lattice-model  (tight-binding)  approximation,  and 
should  be  able  to  extend  the  theoretical  predictions  for  these  systems  beyond  the  current 
limitations.  Figure  3  illustrates  the  results  of  this  work. 

(a) 


Figure  3.  (a)  Zero  temperature  ground-state  phase  diagram  for  the  homogeneous  attractive  spin-1/2  Fermi  gas  in 
ID.  The  phase  boundaries  were  obtained  from  the  Bethe  ansatz.  The  points  X,  Y,  T,  and  D.  along  the  circular 
trajectory  C  denote  the  intersections  with  the  phase  boundaries,  (b)  Densities  of  the  majority  spins  n|  (blue), 
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minority  spins  n  J,  (red),  and  their  difference  nf-nj,  (green)  along  the  trajectory  C  in  (a),  as  obtained  from  variational 
states  having  bond  dimensions  D=8,  16,  20  and  the  exact  Bethe  ansatz  result  (solid  curve). 

Liu's  group  also  carried  out  a  study  of  cold  fermions  in  optical  lattices  with  spin 
imbalance  theoretically.  The  predict  the  existence  of  chiral  superfluidity  with  p-wave  symmetry 
from  an  s-wave  interaction.  Chiral  p-wave  superfluids  are  fascinating  topological  quantum 
states  of  matter  that  have  been  found  in  the  liquid  3He-A  phase  and  arguably  in  the  electronic 
Sr2Ru04  superconductor.  They  are  fundamentally  related  to  the  fractional  5/2  quantum  Hall 
state,  which  supports  fractional  exotic  excitations.  Past  studies  show  that  they  require  spin- 
triplet  pairing  of  fermions  by  p-wave  interaction.  Here  the  key  idea  Liu's  group  proposed  is  that 
pairing  of  fermions  takes  place  between  s  and  p  orbitals  of  the  optical  lattice  (see  Fig.  4),  which 
have  even  and  odd  parity,  respectively.  A  crucial  condition  for  realizing  such  pairing  is  to  create 
sufficiently  large  spin  imbalance  so  that  the  spin  up  and  down  Fermi  levels  reside  in  the  s  and  p 
orbital  bands,  respectively,  with  total  fermion  numbers  N_up  <  N_down.  This  p-wave  state  is 
conceptually  distinct  from  all  previous  conventional  p-wave  states  as  it  is  for  the  center-of-mass 
motion,  instead  of  the  relative  motion.  It  leads  to  spontaneous  generation  of  angular  momentum, 
finite  Chern  numbers  and  topologically  protected  chiral  fermionic  zero  modes  bounded  to 
domain  walls,  all  occurring  at  a  higher  critical  temperature  in  relative  scales,  due  to  the  5-wave 
pairing  nature  in  relative  motion. 


a  b  d 


Conventional  p+ip 


Center  of  mass  p+ip 


Figure  4.  (a)  Schematic  picture  of  a  2D  spin-dependent  optical  lattice,  where  the  spin  up  (s  orbital  band)  and  down 
(p  orbital  band)  component  lying  within  different  geometry  lattice  potential,  respectively.  Here  A  and  B  stand  for 
two  different  sites  in  one  unit  cell,~ex  and~ey  are  the  primitive  unit  vectors,  (b  and  c)  Schematic  views  illustrate 
tunneling  tO,  tl,  t2  and  t3  of  fermions  prepared  in  the  s  and  p  orbitals,  respectively,  (d)  Zero -temperature  phase 
diagram  when  tl/t0!48,  t2/t0142 — The  solid  line  illustrates  the  phase  transition  from  normal  gas  (NG)  to  superfluid 
state  (SF).  When  U/t0o7,  the  critical  value  of  t3/t0  as  shown  by  the  dash  line,  beyond  this  threshold  a  phase 
transition  from  px+py  to  px+ipy  superfluid  state  occurs.  When  U/tOZ7,  px+ipy  superfluid  state  is  the  ground  state 
with  non -zero  t3.  The  thick  solid  line  stands  for  a  two-component  superfluid  state,  (e  and  f)  Pictorially  illustrate  the 
distinction  between  centre-of-mass  and  conventional  p |t i p  pairings.  Here,  the  character  of  pjtip  pairing  states  is 
demonstrated  by  the  phase  of  the  superconducting  gap  labelled  by  the  color  of  the  torus.  In  e,  the  Cooper  pair  is  a 
spin-singlet  pair  composed  of  two  spin  species  fermions  and  has  finite  centre-of-mass  orbital  angular  momentum, 
whereas  in  f  it  is  a  conventional  spin-triplet  pair  formed  by  single  species  fermions  and  has  conventional  relative 
motion  angular  momentum. 
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Dynamics  and  Out-of-Equilibrium  phenomena 

If  one  wants  to  experimentally  study  an  equilibrium  phase  diagram,  the  system  must  be  in 
equilibrium.  During  the  OLE  project,  we  realized  that  this  tautology  represented  one  of  the  most 
challenging  aspects  of  the  OLE  program.  Great  efforts  have  been  made  to  understand  the 
dynamic  processes  which  drive  an  optical  lattice  system  to  equilibrium. 

DeMarco  and  Ceperley  collaborated  on  two  results  involving  dynamics  and  out-of¬ 
equilibrium  phenomena  in  strongly  correlated  optical  lattices.  The  first  focused  on 
measurements  and  simulations  of  the  peak  fraction  for  a  Bose-Einstein  condensate  loaded  into  a 
cubic  optical  lattice.  DeMarco’s  group  measured  how  peak  fraction  varied  as  the  entropy  per 
particle  of  the  gas  and  the  lattice  potential  depth  were  varied  in  the  superfluid  regime  (Fig.  5). 
These  measurements  were  compared  with  large-scale,  exact  quantum  Monte  Carlo  simulations 
carried  out  by  Ceperley’s  group.  A  discrepancy  was  discovered  at  high  temperatures — the 
experimentally  measured  peak  fraction  was  too  high  for  a  gas  in  thermal  equilibrium  in  the 
lattice.  An  explanation  for  this  phenomenon  in  terms  of  suppression  of  Landau  damping  (Fig.  6) 
was  found  in  collaboration  with  Stefan  Natu  at  the  JQI.  These  results  were  published  in  Physical 
Review  A  91,  023625  (2015)  and  highlighted  as  an  Editor’s  Suggestion.  The  suppression  of 
thermalization  observed  in  this  system  has  critical  implications  for  optical  lattice  experiments. 


S/A/  (kB) 


Figure  5.  Images  taken  of  ultracold  atoms  prepared  at  difference  entropies  per  particle  S/N  after  turning 
on  an  optical  lattice  to  depth  s. 
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Figure  6.  Schematic  diagram  of  Landau  damping,  whereby  a  quasiparticle  (y)  can  decay  to  two  particles 
(p  and  v). 

DeMarco’s  and  Ceperley’s  groups  also  collaborated  on  a  project  to  understand  the  analog 
of  the  Kibble-Zurek  scenario  for  a  disorder  quench.  DeMarco’s  group  measured  excitations 
generated  by  quenching  disorder  for  a  superfluid  loaded  into  a  disordered  optical  lattice  (Fig.  7). 
A  threshold  disorder  required  to  generate  excitations  was  found  and  compared  with  the  disorder 
required  to  create  a  Bose-glass  phase  in  exact  quantum  Monte  Carlo  simulations  carried  by 
Ceperley’s  group.  These  QMC  simulations  are  the  most  complex  and  largest  scale  ever  carried 
out  for  a  disordered  system.  Both  qualitative  and  quantitative  agreement  between  the  onset  of 
excitations  in  the  experiment  and  Bose-glass  in  the  simulations  was  demonstrated  (Fig.  8).  This 
behavior  is  evidence  that  the  Kibble-Zurek  mechanism  may  apply  to  quantum  phase  transitions 
in  disordered  systems,  which  is  an  important  issue  for  approaches  such  as  adiabatic  quantum 
computing.  This  work  is  available  as  a  preprint  (arXiv:  1503.02333  (2015)). 


Figure  7.  Quenching  procedure  carried  out  by  DeMarco’s  group.  The  disorder  strength  A  is  rapidly  reduced.  For 
sufficient  disorder,  excitations  such  as  vortices  are  created  during  the  quench  (i),  while  excitations  are  absent  at 
lower  A  (ii). 
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Figure  8.  Comparison  of  excitations  measured  in  the  experiment  (a)  and  the  fraction  of  Bose-glass  atoms  in  the 
simulations  (c).  The  inset  to  (c)  shows  the  lowest  two  occupation  eigenfunctions  of  the  single -particle  density 
matrix  for  A=0  (left)  and  A=1  ER  (right).  Images  of  excitations  generated  in  the  experiment  are  shown  in  (b),  and 
simulations  of  the  occupation  for  atoms  in  superfluid  (blue)  and  Bose-glass  phases  (green)  are  shown  in  (d). 


DeMarco  and  Mueller  collaborated  on  an  effort  to  understand  rethermalization  for  bosons 
trapped  in  a  strongly  correlated  optical  lattice.  They  found  that  phase  coherence  is  more  robust 
than  one  would  naively  expect.  One  can  maintain  equilibrium  during  ms  ramps,  even  though  the 
characteristic  hopping  time  is  an  order  of  magnitude  larger.  This  is  good  news  for  exploring 
equilibrium  physics,  and  is  due  to  the  role  of  interactions.  DeMarco’s  group  developed  a  method 
for  rapidly  bringing  the  quasimomentum  distribution  of  the  gas  out  of  equilibrium  (Fig.  9).  The 
re-equilibration  time  of  a  thermal  gas  was  measured  for  different  lattice  potential  depths  in  the 
superfluid  regime.  Rapid  thermalization  and  a  strong  violation  of  the  Mott-Ioffe-Regel  criterion 
was  observed,  which  rules  out  Boltzmann  transport  approaches  to  understanding  thermalization. 
Similar  effects  have  been  observed,  but  not  understood,  in  strongly  correlated  electron  solids. 
Mueller  modeled  the  rethermalization  process  using  a  simple  quantum  approach.  Agreement 
with  no  free  parameters  between  the  measurements  and  model  was  found  for  high  lattice  depths 
(Fig.  10).  These  measurements  are  important  to  efforts  to  cool  atoms  in  lattices  and  may  shed 
light  on  unusual  diffusion  and  transport  phenomena  in  solids.  These  results  are  available  as  a 
preprint  (arXiv:  1503.07606  (2015))  and  are  under  review  for  publication  in  Nature  Physics. 
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Figure  9.  Procedure  to  disrupt  equilibrium  using  stimulated  Raman  transitions  developed  by  DeMarco’s  group 
(a).  Images  of  the  gas  taken  after  turning  on  the  lattice  (top)  and  after  removing  atoms  (bottom)  are  shown  in  (b). 
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Figure  10.  Measured  rethermalization  rate  x  normalized  to  the  tunneling  time  h/t  for  different  lattice  potential 
depths  s.  The  quantum  model  of  thermalization  is  shown  as  dashed  line. 


Other  Theoretical  Issues 

In  addition  to  those  mentioned  above,  our  team  have  explored  various  theoretical  issues 
that  are  not  necessarily  directly  tied  to  the  experiments  carried  out  in  the  experimental  groups 
within  our  team,  but  are  broadly  related  to  the  OLE  program. 

Several  groups  have  explored  topological  phases  in  ultracold  atoms  subject  to  spin-orbit 
coupling  (SOC).  Pu'  group  has  investigated  properties  of  both  bosonic  and  fermionic  systems 
with  spin-orbit  coupling.  They  also  proposed  to  create  SOC  with  quantized  light  field  supported 
by  an  optical  cavity,  where  the  atomic  feedback  onto  the  cavity  field  give  rises  to  a  dynamic 
SOC.  Mueller's  group  has  investigated  ways  of  studying  "topological  phases"  in  an  optical 
lattice  experiment.  These  are  exotic  states  which  may  have  application  in  spintronics  or  quantum 
computing  —  and  which  have  become  an  area  of  intense  study  over  the  past  decade.  They 
conducted  theoretical  investigations  into  the  possibility  of  observing  these  states  in  cold  atom 
experiments.  At  the  moment  there  are  still  technical  difficulties  which  make  such  experiments 
challenging,  but  these  works  (and  related  ones  from  other  authors)  show  that  cold  atom 
experiments  have  advantages  over  solid-state  experiments  for  directly  observing  some  of  the 
more  exotic  features  of  these  states.  For  example,  they  showed  that  radio  waves  can  manipulate 
the  non-local  quantum  correlations  that  these  states  possess.  Ho's  group  has  also  studied 
interaction  effects  in  spin-orbit  coupled  quantum  gases,  and  larger  spin  fermions  with  SOC. 

They  have  shown  that  spin-orbit  coupled  large  spin  fermions  possess  very  rich  spin  structures, 
some  of  which  have  the  symmetry  of  a  Platonic  solid,  and  many  have  non- Abelian  line  defects. 

In  addition,  they  continue  to  make  significant  headway  in  the  following  problems:  (a)  cooling  of 
the  spin  degrees  of  freedom  of  lattice  quantum  gases,  and  (b)  topological  matter  caused  by 
external  fields  and  by  interaction  effects. 


10 


Mueller's  group  also  carried  out  studies  of  the  feasibility  of  observing  Nagaoka 
ferromagnetism  in  an  optical  lattice  [Yariv  Yanay,  Erich  Mueller,  Veit  Elser,  Magnetic  polarons 
in  two-component  hard  core  bosons,  Phys.  Rev.  A  87,  043622  (2013)  ].  This  is  a  phase  of  an 
iconic  magnetic  model  which  was  used  in  the  60 ’s  to  develop  understanding  about  various 
mechanisms  for  ferromagnetism.  Experimentally  studying  this  phase  is  an  important  piece  of 
validating  theoretical  techniques  which  are  being  applied  to  transition-metal-oxides.  They  found 
that  a  “quantum  gas  microscope”  experiment  on  cold  bosons  could  readily  observe  the 
development  of  ferromagnetic  order. 

Several  groups  have  extended  this  DARPA  funded  research  into  other  directions,  due 
to  the  questions  from  the  leaders  in  other  communities,  Ho  has  begun  to  think  about  extending 
the  Quantum  Simulation  Program  of  cold  atoms  to  other  disciplines  (from  the  traditional  ultra- 
low  temperature  Helium  physics  to  lattice  gauge  field  theory  in  high  energy  physics),  as  well  as 
the  common  challenges  of  these  different  fields.  For  example,  he  has  started  new  investigations 
of  quantum  gases  in  curved  space  and  the  highly  non-equilibrium  phenomena  of  inflationary 
expansion  of  quantum  gases.  In  a  related  study,  Mueller  has  found  that  when  one  changes 
parameters  faster  than  the  equilibration  times,  “coherence  waves”  spread  ballistically  through  the 
atomic  cloud  (see  Fig.  11).  As  pointed  out  by  Hung  et  al.  [Science  341,  1213  (2013)],  these 
waves  are  visible  in  experiments,  and  are  analogous  to  structures  seen  in  spatial  maps  of  the 
black-body  radiation  spectrum:  In  the  expansion  phase  of  the  early  universe,  the  size  of  the 
universe  changed  faster  than  any  equilibration  times. 
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Figure  11.  Illustration  of  coherence  wave  following  a  quench  into  the  superfluid  state. 
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